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Abstract 

o\ ■ 

■ The behaviour of the Maxwell field near one of the spatial infinities of the Schwarzschild 

solution is analysed by means of the transport equations implied by the Maxwell equations 
on the cylinder at spatial infinity. Initial data for the Maxwell equations will be assumed to 
^ ■ be expandable in terms of powers of a coordinate p measuring the geodesic distance to spatial 

t* — ' infinity (in the conformal picture) and such that the highest possible spherical harmonics at 

order p are 2 p -polar ones. It is shown that if the 2 p -polar harmonics at order p in the initial 
data satisfy a certain regularity condition then the solutions to the transport equations at 
orders p and p+1 are completely regular at the critical sets where null infinity touches spatial 
infinity. On the other hand, the solutions to the transport equations of order p + 2 contain, 
in general, logarithmic singularities at the critical sets. It is expected that the ideas and 
techniques developed in the analysis of this problem can be employed to discuss the more 
challenging case of the transport equations implied by the conformal Einstein equations, and 
thus provide a way of constructing a proof of a certain rigidity conjecture concerning the 
5— i , developments of conformally flat initial data sets. 
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1 Introduction 



In |13j the following rigidity conjecture concerning the development of conformally flat initial 
data for the Einstein vacuum equations sets was put forward: 

Conjecture. If an initial data set for the Einstein vacuum equations which is time symmetric 
and conformally fiat in a neighbourhood B a {i) of infinity yields a development with a smooth null 
infinity, then the initial data set is Schwarzschildean in B a (i). 

There is further evidence for generalisations of this conjecture to more general classes of initial 
data — see e.g. [HI Q31 [TS]. A presumptive proof of such a conjecture would entail two aspects: 

(I) A deep understanding of the mechanisms of the asymptotic expansions that are obtained 
using Friedrich's construction of the cylinder at spatial infinity — see e.g. [31 [S] • 

(II) The construction of estimates which allow to get around the manifest degeneracy of the 
Einstein conformal equations at the critical sets were null infinity "touches" spatial infinity 
— see e.g. [J. 

The present article explores a set of ideas on how to deal with the aspect (I) mentioned above. 

Friedrich's seminal work on the i° problem has rendered a representation of spatial infinity 
which exploits to the maximum the conformal properties of the spacetime by choosing a convenient 
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gauge. In particular, spatial infinity is represented in this gauge by an extended set, the cylinder 
at spatial infinity, X — see [3] . This framework uses the extended conformal field equations — as 
given in, say, [3J [3J — to pose a regular initial value problem at spatial infinity. 

The cylinder at spatial infinity can be regarded as a limit set of incoming and outgoing 
light cones. It has the remarkable property of being a total characteristic of the conformal 
Einstein propagation equations. This property implies that the propagation equations reduce in 
their entirety to an interior system of transport equations upon evaluation at X. This allows to 
transport information from some initial data prescribed on some initial hypersurface S up to the 
critical sets I where spatial infinity touches null infinity. Let = (0o, 01 j 02, 03, 04) denote the 
components of the rescaled Weyl tensor in a particular gauge and with respect to a particular 
frame — to be described in the subsequent sections. Then the transport equations at X allow to 
calculate d^<p\x — where p is some radial coordinate measuring the distance to X — on the whole 
of X — and in particular at X ± — if the value of <9^0 onl = X (IS is given and the lower order 
derivatives dp<f)\x, k — 0, ... ,p — 1 are known. This unfolding of the evolution process can be 
thought of as the construction of the Taylor like expansions of the form 

p 

A general observation is that the derivatives dP<j)\x contain logarithmic divergences at X ± . By 
logarithmic divergences at, say r = ±1, it will be understood the appearance expressions of the 
form 

(l-T)Mn(l-r), (l + r) fe +ln(l + r). (2) 

Note how the exponents fc_ and k + regulate the smoothness of the term. The logarithmic diver- 
gences appearing in the solutions of the transport equations implied by the conformal Einstein 
equations can be grouped in two types: 

(i) those appearing from general linear aspects of the equations; 

(ii) and those arising from the particular nonlinear aspects of the equations. 

The logarithmic divergences of (i) are well understood and their preclusion imposes mild 
conditions on the initial data — see [3J [TT]. Not surprisingly, it is the understanding of the 
logarithmic divergences of type (ii) which poses the real challenge. The discovery of a particular 
class of this logarithms led to conjectureQ] At the moment, there are no general results regarding 
the behaviour of this kind of logarithms. The only knowledge available is what has been inferred 
from the calculation of the expansions ([I]) up to a finite order — say p — 7, 8 — using computer 
algebra methods. 

With regards to the coefficients (<9^0|z) in the expansions ([1]), one would like to prove a result 
on the lines of: 

Subconjecture. Given a time symmetric initial data which is conformally flat on B a {i), and if 
dp<P\i ar & smooth at X for allp, then 3 p 0|x° is Schwazschildean. 

In other words, if one starts with a general time symmetric, conformally flat initial data set 
and one works out the behaviour of the derivatives <9p0|x at X , then one will observe that 
generically there are logarithmic terms. If one removes the parts of the initial data giving rise 
to these divergences at every order p, then one will observe that the jet {9 p 0|xo, p = 0, 1, . . .} 
is equal to the one obtained from looking at the t = constant slice of the Schwarzschild solution 
in isotropic coordinates. In order to make these ideas work, one needs the understanding of the 
asymptotic expansions which was alluded to in point (I). 

The complexity and vast number of equations which has to be controlled when dealing with the 
conformal Einstein equations makes developing the required understanding a daunting endeavour. 
That is the reason why it is desirable to have simpler "toy models" on which diverse techniques 
could be tested and from which underlying structures could be apprehended. Devising such a toy 
model is a delicate issue. Of course, one wishes for a tractable model, but at the same time one 
would like to avoid simplifying things too much and render the analysis trivial. 
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The obvious "toy model" in the form of a spin-2 zero-rest-mass field propagating close to one of 
the spatial infinities of the Schwarzschild spacetime is not an adequate one in view of the existence 
of an obstruction in the form of the Buchdahl constraint which precludes the consideration of 
generic fields in the analysis — see e.g. pQ. This problem is overdetermined. 

Due to the above considerations, attention is turned to another possible toy model: the 
propagation of the Maxwell field (spin-1 zero- rest- mass field) in the vicinity of one of the spatial 
infinities of the Schwarzschild spacetime. A direct application of the ideas of, say, [11] shows the 
existence of logarithmic divergences at the sets T^- of the underlying Schwarzschild spacetime 
which are due to generic linear features of the Maxwell equations which are shared with the 
conformal Einstein equations — the type (i) divergences. An understanding of the mechanism 
behind these divergences, as it will be seen, poses no challenge. The real difficulties arise from 
the existence of another type of logarithmic divergences which are generated by the interaction 
of the Maxwell field with the background spacetime. A preliminary analysis of the analogue of 
the expansions ^ for this toy model using the computer algebra methods of [T3] reveals that the 
inner workings of the logarithmic divergences are similar to those of the nonlinear propagation of 
the gravitational field. 

The main goal of the present work is to develop a procedure of tracking down which parts 
of the initial data are responsible for specific logarithmic divergences at the critical sets. To this 
end it shall be assumed that one has a class of initial data for the Maxwell equations expandable 
in powers of a geodesic distance p to infinity. Further, it will be assumed that the coefficients of 
this expansion at say, order p p admit at most 2 p -polar harmonics — this happens if, for example, 
the initial data is analytic in a neighbourhood of infinity. Consistent with the discussion of the 
previous paragraphs, it will also be assumed that the initial data satisfies the regularity conditions 
precluding the appearance of logarithms of the so-called type (i). Given this class of initial data, 
the question is the following: at which order — if at all — do the 2 p -polar terms in the initial data 
produce divergences at the critical sets? Transport equations for 2 p -polar quantities appear for 
the first time at order p. The present analysis shows that the solutions at orders p and p+ 1 do not 
contain logarithmic divergences — for order p this is direct, but for order p+ 1 there is an intriguing 
"conspiracy" of terms taking place. Moreover, for order p + 2 the solutions contain, in general, 
the logarithmic terms. A similar pattern — although with a further layer of complication — arises 
for the conformal Einstein field equations. Hence the interest of the present analysis. 

It should be mentioned that there has recently been some related work on the peeling proper- 
ties of solutions of the wave equation on a Schwarzschild background — see [8] . These analyses 
can, most certainly, be extended to the case of the Maxwell equations. In any case, the scope, 
aims, techniques and motivations for these results are different to the ones in this article. 

The article is structured as follows: section 2 briefly discusses briefly the Schwarzschild space- 
time in the so-called F-gauge in which the construction of the cylinder at spatial infinity is carried 
out. Section 3 is concerned with the Maxwell equations in the present context. Section 4 is consid- 
ers further properties of the transport equations implied by the Maxwell equations at the cylinder 
at spatial infinity, including some regularity conditions. Sections 5, 6 and 7 discuss, respectively, 
the solutions to the order p, p + 1 and p + 2 transport equations. There are some concluding 
remarks in section 8. Finally, there are three appendices including the definitions of some spino- 
rial objects, the first orders of the asymptotic expansions of the Schwarzschild solution, and some 
properties of Jacobi polynomials which are used in the analysis. 

2 The Schwarzschild spacetime in the F-gauge 

The Schwarzschild spacetime was discussed in reference [3] as an example of the implementation 
of the framework of the cylinder at spatial infinity — or F-gauge. This representation will be used 
to analyse the propagation of the Maxwell field near spatial infinity. Some of the relevant features 
of this construction are reviewed. 

The Schwarzschild metric in isotropic coordinates is given by 
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Now, consider the time symmetric hypersurface S. Writing h = Q h, one can introduce a 
conformal metric and conformal factor given by 

h = -W + f6 *) i V = J^JW> (4) 

where the radial coordinate p — 1/r has been introduced. Let i\ and i% denote the infinities 
corresponding to the two asymptotic ends of the hypersurface S. Further, let S = S U {«i,«2}- 
Let i — i\ denote the infinity corresponding to the locus p = 0. The discussion in this article will 
be concerned with the domain of influence, J + (B a (i)) of a sufficiently small ball, B a (i), of radius 
a based on i. The point i can be blown up to a 2-sphere, S 2 . With this idea in mind introduce 
the set C a = {B a (i)\i)US 2 . 

The use of a gauge based on conformal Gaussian coordinates to discuss the structure of B a {i) 
leads to a conformal factor for the portion of spacetime under consideration which is given by 

e = -(i-r 2 4), (5) 

where 

u= , = = p(l + pm/2), (6) 

and k > is a smooth function such that n = kp with k{i) — 1. The function k encodes the 
remaining conformal freedom in the setting. In what follows, and in order to ease the calculations 
let 

k = p. (7) 

The coordinate r is an affine parameter of conformal geodesies whose tangent at r = is parallel 
to the normal of S. The coordinate p can be extended off S by requiring it to be constant along 
the aforementioned conformal geodesies. "Angular coordinates" can be extended in a similar 
fashion. 

In view of the conformal factor (J5]) define the manifold 

M a . K = {(T,q) \qeC a , -- <T< -}, (8) 
such that J + (B a (i)) C M. a ,Ki an d the following relevant subsets thereof: 

1 = {(r, g) e -Ma.K | P (q) = 0, |r| < ^} , (9a) 

X ± = {(r, g) S M a , K | p(q) = 0, r = ±1} , (9b) 
l° = {{T,q)€M a , K \p = 0, r = 0}, (9c) 

S ± = {(T,q)e M a<K \qeB a (i)\I°, ^ = ±^}, (9d) 

denoting, respectively, the cylinder at spatial infinity, the critical sets where spatial infinity 
touches null infinity, the intersection of the cylinder at spatial infinity with the initial hyper- 
surface S, and the two components of null infinity. 

It will be convenient to work with a space-spinor formalism — see [9] . In order to write down 
the field equations, introduce a null frame caa' satisfying g(cAA>,CBB') = e AB^A'B'- Let t A a* 
— with normalisation t aa t^a' — 2 — be tangent to the conformal geodesies of which r is a 
parameter. The frame can be split in the form 

CAA' = ^jT AA ,T CC C C C ~ T B A ,C AB , (10) 

with 

t AA c AA > = \/2d T , c AB = T^ A B c B)B ,. (11) 
For the purpose of the present investigation, the following choice is made: 

coo' = ±={{l + c Q )d T + c x d p ), c iv = -J= ((1 - c°)d T - Jdp) . (12) 
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The remaining vectors of the frame, coi' and c 1( y must then be tangent to the spheres {t = 
constant, p — constant}, and thus cannot define smooth vector fields everywhere. To avoid 
this difficulty all possible tangent vectors Cov and cic will be considered. This results in a 5- 
dimensional submanifold of the bundle of normalised spin frames. Rotations coy — > e lQ , a G M 
leave this submanifold invariant. Hence it defines a subbundle with structure group U(l) which 
projects into M UyK . All the relevant structures will be lifted to the subbundle, which in an abuse 
of notation will be again denoted by M a , K - 

The introduction of the bundle space M a . K and of the frame caa> in our formalism implies 
that all relevant quantities have a definite spin weight and hence admit a characteristic — spherical 
harmonics — expansion in terms of some functions T- k l associated with unitary representations 
of SU(2, C) — see e.g. [3j [5] for a more detailed discussion in this respect. Further, one can 
introduce differential operators X, X + and X- defined by their action on the functions T- { . It 
is noted that with the help of these operators one can write 

caa> = c AA ,df, = c° AA ,d T + c\ A ,d p + c\ A ,X+ + c AA ,X_. (13) 

In addition to the frame caa 1 , the geometry of M a . K is described by means of the associated 
connection Taa'bc, the spinorial counterpart of the Ricci tensor of a Weyl connection Oaa'bb' 
and the rescaled Weyl spinor 4>abcd- The unprimed (space-spinor) version of the connection 
spinor is given by T A bcd — t b b T A b'CD, ^aa'cd = ^abcdt b A '-, which is decomposed as 



FaBCD = ~J= (S.ABCD - X(AB)Cd) - T^AbIcD- (14) 

Similarly, one considers Qabcd = t c a t d b &aa'bb'- The explicit spherical symmetry of the 
spacetime justifies the following Ansatz in terms of irreducible spinors: 

C AB = C °XAB: c\ B = C XABi C AB = C VAB, c\ B = C + Z A B, (15a) 

fAB = fXAB, iABCD = £(tACXBD + i-BD^Ac), (15b) 

X(AB)CD = X2t 2 ABCD + Xhh>ABCD, (15c) 

Qabcd = Qi^abcd + QhhABCD + QxCabxcd, 4>abcd = 4> e ABCD- (15d) 



The definitions of the irreducible spinors introduced above are given in the appendix. The man- 
ifest spherical symmetry of this representation implies that the functions c°, c 1 , c^, f, £, \2, Xh, 
4> have spin- weight 0. Furthermore they only contain the function T ° = 1. 

The functions c°, c 1 , c^, /, ^, X2, Xh, 4> are determined by solving the conformal propagation 
equations discussed in [3] with the appropriate initial data. The problem of reconstructing the 
conformal Schwarzschild solution from the given data amounts to finding a solution u = u(r, p; m) 
of an initial value problem of the type 

d p u = F(u,T,p;m), u(0, p; m) = u a (p; to), (16) 

with analytic functions F and uq. The solution with to = corresponds to a portion of the 
conformal Minkowski spacetime, in which the only non- vanishing components of the solution are 
given by 

c° = -r, c x =p, c± = l, f=l. (17) 

Since in this case the solution exists for all t, p £ R, it follows from standard results of ordinary 
differential equations that for a given m there is sufficiently small po such that there is an analytic 
solution to the system P^|) which extends beyond J' for \p\ < \p \. Hence, one can recover the 
portion of the Schwarzschild spacetime which lies near null and spatial infinity, J + (B a (i)), if a is 
taken to be small enough. 

It follows from the above discussion that the coefficients that are obtained from solving the 
transport propagation equations on the cylinder at spatial infinity — as discussed in [3l fT3]— 
correspond to the first terms in the expansions of the solutions of (|16p . See the appendix for a 
list of the relevant expansions. 
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3 The Maxwell field 



As it is common when working with spinors, consider the Maxwell tensor written as 

FaA'BB' = <t>ABtA'B> + <t>A'B>£AB- (18) 

Then the Maxwell equations are equivalent to 

V AA '<f> AB =0, (19) 

that is, the spin-1 zero-rest-mass field equations. Recall, as well, that if the conformal weight of 
4>ab is chosen properly, the vacuum Maxwell equations are conformally invariant. The spinor 
4>ab is symmetric and thus, one can write 

<i>AB = 0oe% + 01 e AB + H<?ab- (2°) 

In the F-gauge, and using a space-spinor decomposition, the above equation is equivalent to 
the following propagation equations: 

(V2 - 2c° 01 )dr<bo + 2c° o 3 T 0i - 2c^d a <j)o + 2c% d a fa 

= (2r 00 ii - 4rioio)0o + 4riooo0i - 2r OO oo02, (21a) 

V28 T (f>l - C^dr^f) + Co O <9 T 02 ~ C^da^Q + CQ d a (f> 2 

= -(2r mo + fn)4>o + 2(r 110 o + r OO ii)0i - (2r 00 oi - /oo)<fe, (21b) 

{V2 + 2c° Q1 )drh - 2c\ x d T 4>i + 2c^d a dp 2 - 2c? 1 a <fc 

= -2Yi\i\4>o + 4r O in0i + (2rnoo - 4r ioi)</>2- (21c) 

In addition, one has the following constraint equation 

Cnd T (j) - 2cQ 1 d T (f> 1 + CQ Q d T (f> 2 + c^dafo - 2cQ 1 d a cj) 1 + CQ d a (j)2 

= (2r mo — 2r( i)ii)(/'o + 2(r on — rn O o)0i + (2r( i)oo — roooi)^2- (22) 

Because of the symmetries of the Schwarzschild spacetime, the connection coefficients Toon, Fioio, 
Tnoo and Toioi are the only non-vanishing ones. The leading terms of frame and connection 
components are listed in the appendix [Bl 

Writing </> = (<j)Q, 4>i, ^2), the Maxwell propagation equations can be concisely expressed in the 
form 

V2Ed T + A AB c AB d^ = B{Y ABCD )<t>, (23) 

where E denotes the (3 x 3)-unit matrix, A AB c AB are (3 x 3)-matrices depending on the co- 
ordinates and B(T) is a matrix valued function of the connection Tabcd- For consistency the 
constraint equation (|22|) will be written as 

F AB c!\ B d^ = H{Y ABCD )^ (24) 

where F AB c AB denotes (1x3) matrices, and H is a function of the connection. 

The system (|23|) shares the property with the Bianchi propagation system for the rescaled 
Weyl tensor and the propagation system for linearised gravity that its principal degenerates at 
critical sets . More precisely, one has that 

{V2E + A AB c° AB )\ T = \/2diag(l - r, 1, 1 + r). (25) 

Most non-smooth and/or polyhomogeneous behaviour of the solutions of the Maxwell equations 
at the conformal boundary can be traced back to precisely this degeneracy. 

Another deciding property of the Maxwell equations in the aforediscussed form stems from 
the fact that A AB c AB \j — 0. From here it follows — see e.g. the discussion in JT] — that the 
equations imply an intrinsic system of transport equations on X which determines the solution 
on T in terms of the data implied by <\> onl". In other words, X is a total characteristic of (|23p 
and (1211). 
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As in the case of the conformal Einstein field equations of [3], given an asymptotically flat 
spacetime one can produce, for the Maxwell field, a hierarchy of transport equations at the 
cylinder at spatial infinity out of the equations (j2"3"|) and (|24[) . Differentiating the latter equations 
with respect to p, say, p times and evaluating them at the cylinder at spatial infinity T one 
obtains, by virtue of c\ B \x = the equations 

(V2E + A AB (c° AB )^)d r ^ + A AB (c%)^d c ^ 

= B{Tf BCD )^ + £ ( P ) (B(T ABCD )^ A AB (c>X B )^d^y), (26a) 

3=1 ^ 

F AB (c° AB )^d T ^ + F AB {c c AB )^d c ^ 

= H{Tf BCD )^ + £ ( P ) mTABCD)^ - F AB {cf AB )^d^-% (26b) 

3=1 

with C = ±. It is noted that the principal part of the above equations is universal, in the sense 
that it depends only on the zeroth-order solution of the background metric. The first term on the 
right-hand side is also universal, while the remaining non-homogeneous terms depend on lower 
order terms of the Maxwell field and higher order terms of the background spacetime. As in 
the case of the gravitational field, provided some suitable initial data one can solve the above 
equations to any desired order in a recursive way. 

The coefficients 4>o, </>i and 4>2 have, respectively, spin- weights 1, and —I. Consistently with 
their spin-weight, the following Ansatz will be made: 

oo p 2q 

fo= E E H-\ a 3^ k ( T ) T 2 q k q -l +1 P P i ( 2? ) 

P=|i-i|9=|i-j| k =o 

with j = 0, 1, 2, and a^ p -2 q ,k '■ K — > C. 

The above assumption will allow to solve the hierarchy of transport equations (|26a[) and (|26b[) 
by "decomposing in modes" . It is important to note that because the background terms in the 
transport equations only contain the function T ° = 1, it will not be necessary to linearise 
products of the form Tj 1 mi x T- 2 m ^ when performing the mode decomposition. 

3.1 Initial data for the Maxwell equations 

The equation (|22| evaluated on r = reduces, for the Schwarzschild spacetime, to 

pd p ^ + \x^ a - X -X+<\> 2 = 0. (28) 

The latter equation can be read as an equation (ODE) for 0i, if cpo and <j>2 are given on the initial 
hypersurface. In the present work, initial data for the Maxwell equations on the initial hyper- 
surface S which when lifted to C a is expandable in powers of p will be considered. Accordingly, 
prescribe on C a components 4>o and (f>2 of the form 

oo p 2q oo p 2q 

= EEE/»a'(°) t ,V/- = EEE -^Ao)T 2q v> p , (29) 

p=0 q=l fc=0 p=0 ?=1 fc=0 V' 

which correspond to the expansions of analytic <^>o and 4>2 — see [3J. One can use equation (|22[) 
to obtain the coefficients ai, P ;2g,fe(0) in 

oo p 29 

^^EEE/^^V- ( 30 ) 

p=0 9 =0 fc=0 ^' 

in terms of ao,p;29,fc(0) and a2, p; 29,fc(0). It is noted that in this construction the coefficient 
ai,0;0,o(0) — the electric charge — is not determined and can be prescribed freely. On the other 
hand, from equation (f2"2"| it follows that ai, p; o,o(0) = for p > 1. 
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4 Further properties of the transport equations 



The discussion in this section is inspired from the treatment given in [3] of the transport equations 
for the Conformal Einstein equations. Besides being given as a reference, it will help to highlight 
the difficulties to be confronted when analysing the solutions of the equations (|26a[) and (|26b[) for 
an arbitrary order p. 

4.1 Reduction to an overdetermined system of ODE's 

The propagation transport equations implied by the Maxwell equations on the cylinder at spatial 
infinity of the Schwarzschild background are of the form: 

(1 + r)d T ^ + X + 4 P) - (p - 1)^ = i#>, (31a) 

d T ^ ] + \{X + ^ + X^ ] ) = B$\ (31b) 

(1 - r)d T 4 p) + + (p - 1)4 P) = 4 P \ (31c) 

while the constraint transport equation is given by 

rdj? + \x + <t>¥ \x^ = S (P) - (32) 

As mentioned before, the terms R^ , R^ and S 1 ^ depend solely on the lowers order 

solutions to the equations. Using the expansion Ansatz (j27|) it can be seen that the coefficients 
&j,p;q,k satisfy the equations 



(1 + T)a' 0iP . q>k + yjq{q + l)ai^ q . k - (p - l)a 0)P - q ,k = U 0)P;q}k , (33a) 
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a 'l,p;q,k + ^V^O? + l)(02,p;g,fc ~ ao,p;q,k) = Ul, p -q t k, (33b) 



(1 - T ) a 2,p;q : k - V<l{Q + l)ai,p;q,k + (P~ l)a2,p;q,k = ^2,p;g,fc, (33c) 

and 

Ta hp;q,k + 2^9(9+ l)(02,p;g,ft + a 0,p:q,k) ~ P0Ll,p;q,k = T p ;q.k, (34) 

for p > 1, 1 < q < p, and where again the nonhomogeneous terms f/o,p;g,fe> Ui, p - q .k, U 2 ,p-,q,k 
and T p - qi k in the equations depend on the background spacetime and lower order (Jj,p;g,fe's. The 
problem has been reduced to the analysis of an overdetermined system of ordinary differential 
equations. The last equation can be used to eliminate the unknown ai lP;9j fc. In what remains of 
this section, attention will be focused on a fixed — but otherwise arbitrary — value of the indices 
p, q, k in their range of definition which henceforth will often be suppressed in order to avoid 
cluttering the equations. 

The coefficient a\ can be eliminated from the following considerations by noting that 

1 r~, w , 1 



- -V<?(<?+1)(1 + r)a - ~V<?(g+l)(l - r)a 2 + pa x = tU x - T, (35) 
from where one can solve for a±. One is left with the system 

(1 + T)a' Q + + 1)(1 + r) - p + lj a + ^-q(q + 1)(1 - r)a 2 



= U + - ^q{q + l)(T -tU x ), (36a) 



(l-T)4-i-g(g+l)(l + r)a + (~q(q + 1)0- - r) +p - 1 j 



= C7 2 + -V9(9+1)(tC/i-T). (36b) 



For future reference the above system will be written in matricial form as 

Y p.iA t ) = A p , q .k( T )Y p , q . k {T) + B p , q . k {T), 



(37) 



with Ap. iQt k a (2 x 2)-matrix, and Y p - q ^, B Ptq k (2 x l)-column vectors. Note that the Wronskian 
determinant detX p;g of the above system — satisfying (detX p . q )' = (TrA p . qt k)X p - q — is given by 

detX m =X! p ) . q (l-T 2 y- 1 , (38) 

where X® is a constant. In order to calculate the fundamental matrix of the system consider 
the equations (again, the multisubindex begin supressed) 



(1 - r 2 K + (2 - (2 - 2p)r)a' + {q +p)(q - p + 1) = 0, 
(1 - r 2 K + ( - 2 - (2 - 2p)r)a' 2 + (q + p)(q - p + 1) = 0, 



(39a) 
(39b) 



which arc implied by the homogeneous parts of the equations (|36ap and (|36b[) . Equations (|39a[) 
and (|39b[) arc so-called Jacobi equations 



D( n , a ,0)a = (1 - r)a" + ((3-a~(a + j3 + 2)t)o! + n(n + a + j3 + l)a = 0, 



(40) 



with 



a = — 1— p, f3 = l— p, ni=p + q, n 2 =p-q-l. (41) 

Note that if Oo(t) is a solution of equation (|39a|) . then q,q(t) = a (— r) is a solution of equation 
(|39bj) . A class of solutions is given by the Jacobi polynomials P^^\t) — see the appendix for 
their definition and relevant properties. 

Using the definition of the polynomials, P p+q ~ P ' P are seen to vanish identically. On the 



other hand, 



Q2 - t^p+q 



(42) 



is a solution of degree n 2 = p — q — 1 for p > 2, 1 < q < p — 1. Using some identities of the 
polynomials one finds that 



Qi = 



1 - T 



p+1 



(p+1, -p+1) 



PUT 



(43) 



is another solution, p > 2, 1 < q < p— 1. Clearly Q\ and Q 2 arc linearly independent as n± > n 2 . 
Another solution is given by 



Q 3 = 



l+T 



R 



(-p-i.p-i) 



9+1 



(44) 



for p>2, l<q<p — 1. A direct calculation shows that the fundamental matrix, X Ptq of the 
system is given by 

Qi (-i) g+1 Q 3 N 



x 



(45) 



for p > 2, 1 < q < p — 1. The notation Q{{t) = Qi(— t), Q|(t) = Qz{— r) has been used again. 
Note that the fundamental matrix does not depend on k. The inverse of X p>q is given by 



X 



1 



P,9 y0 

p;<? 



/ (1 + r) 2 p (p+i,i- P ), x (-l) g p (-p-i,p-i) r x 
(l-r)P- 1 9- 1 ^ J (I-t)^ 1 9+1 V) 
(~l) g p (-p-i,p-i) 



'HO 



\(l + r)^ T 9+1 v ■' (1 + r) 

The case p = q, p > 1 is solved separately. It is noted that 

p-i 



0- ~ T ) p(p+M-p) 



p^T^q-l 



(r) 



(46) 



(2p,-p-l,-p+l) 



i_ r V +1 /i + r 



ftfl 



(0,p+l,p-l) a 0j 



(47) 
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where 

£>(o, P+ i, P -i)a = (1 - t 2 )4 + (-2 - (2p + 2))a' , (48) 

so that 

'1 \ P+l /i i \ P— 1 / /-r j 

1 — r \ /1 + T \ / ^ ^ / ds 



— I— J I— J [ C ' +C 'L rai^ j' (49) 

with Ci and C2 some integration constants. Noting, again, that 02 (r) = cto(r) = a (— r) one 
obtains the following 

Proposition 1. For p = q, the fundamental matrix X p . q of the system &36a\) and \36b\) is 
polynomial if and only if 

aO,p;q,k(0) = a2,p;q,fc(0). (50) 

The latter is the analogue, in our context, of the regularity conditions of [3] and [TTj . 

In what follows, it will be assumed that the initial data for the propagation equations i21a\) - 
h21c\) is constructed such that condition \5b)) holds. 

4.2 The solution to the nonhomogenous part 

Once the fundamental matrix of the system (|36a[) - (|36b[) is known, the solution is given by 



a0,p;q,k(l~) 
0>2,p;q,k(T) 



X v<q {r) 



1(0) /«o, P;9 , fc (o)\ + r x -i {s)Bpq{s)ds 

V a 2,p; 9: fc(0)y Jo 



(51) 



As seen before, the entries of X p . q are rational functions with — in the worst of cases — poles of 
order p — 1 at t = ±1. One can make the following direct but crucial statement: 

Proposition 2. If the functions Uo tP - qj k, Ui,p;q,ki U2, p - q ,k and T p - q ^ in equations H21a\) - [^Tc]) 
and A22\) are polynomial in r — that is, if the lower order functions flj,p';g',fe'( T ) with < p' < p, 
> q' > p' , > k' > 2q' are polynomials in r — then the solutions given by &51\) consist of 
polynomials in r and of products of polynomials in r times the functions ln(l — r) and ln(l + r). 

Proof. The entries of the vector X~ q B p>q are, respectively, of the form 

fo(r) , Mj) 

and n Z\7i I » ( 52 ) 



(l + r)(l -t)p~ 1 (l-r)(l + r)P- 1! 

where /o(t) and fzfr) are polynomials in r. Decomposition in partial fractions of the first entry 
renders an expression containing (1 + t) -1 , (1 — t) _ ™ with n = 1, . . . ,p — 1, and a polynomial. 
Analogously for the second entry, where one finds a decomposition in terms of ( 1 — t) ~ 1 , (1 + t)~™ 
with n = 1, . . . , p— 1, and a polynomial. The integration produces, in the first entry, an expression 
containing ln(l ±r), (1 — r) 1 ^™ with n = 2, ...,|>— 1 and a polynomial. Similarly for the second 
entry. Now, the multiplication by the fundamental matrix cancels out the rational expressions. 



Remark 1. From the previous discussion it is clear that even if the regularity condition 
(f50|) holds, the solutions to the transport equations will contain ln(l ± r) terms unless some 
cancellations take place. 

Remark 2. On the other hand, one could use an argument analogous to that of the above 
proof to show that if the lower order coefficients already contain ln(l ± r)-singular terms — say 
because the regularity condition (|50|) is not satisfied, then the integration in equation (|5ip would 
render ln 2 (l ± r)-terms. This type of analysis is, however, beyond the present discussion. 
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5 Expansions at order p 



From the hyperbolicity and linearity of the propagation equations, if for a given order p a certain 
type of harmonics has not yet appeared in the initial data, it will also not be present in the 
solutions to the transport equations. For the class of initial data under consideration, terms 
containing harmonics of the form T 2 1+ . — which will be referred to as 2 p -polar — appear 
firstly at order p. It follows that if p = q, then 

B P ; P Mr) = 0, (53) 

for the class of data under consideration — for all the lower terms ajy-y^' with < p' < p, 
|1 — j\ < l' < p' an( l < k' < 2q' can only contain at most T 2 ^ p _^ k p _ 2+ j functions. In the sequel, 
the following observation will be repeatedly used: the Schwarzschild solution contains only T ° 
sectors. Hence, all the products of coefficients of the Maxwell field and background terms are 
trivially expanded in terms of the functions T 2 2+J - present in the coefficients of the Maxwell 
field. This observation might seem trivial, but note that if one considers a representation of the 
Schwarzschild solution where the spherical symmetry is not manifest, a "mixing of modes" will 
occur — and in particular higher order harmonics would be generated, severely complicating the 
analysis. From here it follows that the first appearance in the solutions to the transport equations 
of terms containing T 2p 1+ . renders no logarithms if the regularity condition (|50[) is satisfied. 

Proposition 3. If the regularity condition H50\) holds, then the solutions ao, P - P ,k(T), ai, p - tP ,k(T) 
and a2, P ; P ,k(T~) of the transport equations h33a\) - ^33c]) and {3$ for p — q and k = 0, . . . , 2p are 
polynomial in r. Moreover, one has that 

a , P ; P ,fc(r) = C p;p>k (l - Tf +1 (l + rf" 1 , (54a) 

ai, p;p , k (r) = C p , p , k ^^(l - T f(l + rf , (54b) 

a 2 , P;P ,fe(r) = C p;Plfe (l - t)p- 1 (1 + t)p+ 1 , (54c) 

with C p - Py k = ao, P ;p,fe(0). 



6 Expansions at order p + 1 



If the condition (|50|) holds, then the first occasion where the harmonics T 2p can produce 

logarithms via the formula (|51[) is in the order p+ 1 of the transport equations. Quite remarkably, 
as it will be seen, this is not the case: cancellations occur which prevent the appearance of 
logarithms at this order. It is necessary to go up to order p + 2 to find the first logarithms in the 
sectors T 2p k p _ 1+ -. More precisely, one has the following: 

Proposition 4. // the regularity condition i5ty) is satisfied for p e N, p > 1 then the functions 
%,p;2p,fe(T) ore polynomial in r — that is, they are logarithm free. 

Remark. It must be emphasised that in view of the formula (|5T|) . the above result arises as 
the result of a remarkable collusion of terms which, in principle, there is no reason to expect. It 
would be very interesting to understand the deeper (possibly group theoretical) reasons behind 
this result. 

Proof. The proof proceeds by an explicit calculation. The equations (I36al) (|36b[) imply upon 
the replacement p — > p + 1, q —> p the system 

, p (1 — r) p (1 — r) , , 

a o, P +i _ 2 (1 + T ) a °' ?)+1 + 2 (1 + T ) a2 ' p+1 ~ °' p+1 ' ^ > 

P (1 + t) p (1 + t) 

a 2, P +i + 2 (i _ r ) a2 -P +1 ~ 2 (1 - T ) a °' p+1 — ^.P+i' \ bbt} ) 
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where the subindex string ;Pi fc has been suppressed throughout for the ease of the presentation. 
In accordance with equation (|26a|) one has 



1 



5 0,p+l 



5 2,p+l 



1+T 
1 

1 - T 



{ho,pO-o,p + h>i !P ai lP + /i2 : pfl2,p + ho,pa'o. p + hi, p a! lp + h2 tP a 2iP ) , (56a) 
(^2, P a o,p + h s l p ai :P + h s p a2, p + h 2p a' 0p + h s lp a' lp + h s 0p a' 2p ) , (56b) 



with hjj, and hj. p for j = 0,1,2 some polynomials in r arising from the order 0,1,2 solutions 
of the transport equations for the Schwarzschild background — see appendix [Bj As before the 
notation W (r) = hj jP (— r) has been used. 

The order p Maxwell transport equations can be used to eliminate a' p , a' lp and a! 2p from 
the expressions for &o,p+i and &2,p+i- Further, equation (|35[) with q = p can be used to eliminate 
a l p+1 . At the end one is left with 



bo.p+i 
i>2,p+i 



(1 + r) 2 
1 

(i~7p 



(/oOO,p + /202 : p) , 
(/l a 0,p + /o°2,p), 



(57a) 
(57b) 



with 



/o(r) 



12 



m(p+l)T((-4p-5)T 6 + (- 



2)r 5 + (25 + 16p)T 4 + (56p + 40)r 3 



-(30 + 48p)7 



24r , 



/ 2 (r) = — p(p + l)m((-4p - 5)r 6 + (4p + 2)t 5 + (29 + 28p)r 4 + (-16p - 8)t 3 



-(18m + 24p)r 2 



It follows then that 



^P+l,P^P:P+l 



Cp;p 
P+l,p 



/a 



\p-i 



-__^fF (l - r) 2 + (-lyjsF; 
H^(l + r) 2 + (-!)*> JJF ) 



(58a) 
(58b) 

(59) 



V (l+r) J 

where the polynomials Ji and J 3 arising from the fundamental matrix are given by 



Jx(r) 



1 



P. 



2p+2 p-i 



(p+2,-p) 



(r), 



1 



Mr) = ^Pp+i 



(-P-2.P) 



(r). 



The polynomials -Fb and F% are given by 

F = (l-r) 2 / + (l + r) 2 /2, 
^2 = (l-r) 2 /| + (l + r) 2 /J, 



(60a) 
(60b) 



(61a) 
(61b) 



In order to determine whether the solutions of the transport equations at this order contain 
logarithms, one has to decide, say, whether 



H(r) = (1 + rf- 1 ((l - t) 2 J-IF + J 3 F 2 ) 



(62) 



and its derivatives have zeros at r = 1. Direct evaluation of formula (|62p using the identities 
listed in the appendix for the values of the Jacobi polynomials at the terminal points r — ±1 
render 



H(l) = 2(1 ~p 2 )m, 

H'(l) = -(p+l)(5m-7pm), 

H"(l) = 0. 



(63a) 
(63b) 
(63c) 
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Using the division algorithm one can write 

H(t) = (1 - T) 2 h 2 (r) + (1 - r)/ii(l) + H(l). (64) 
where /ii(r) and ^(r) are some polynomials. Thus, /12 (1) = 0, if H"(\) = 0, whence 

H(t) = (1 - r) 3 ^^) + (1 - r)/ii(l) + H(l). (65) 

Hence 

C ■ I 3( > ' (1 _ T ) 2 + fi - T ) 3 I 

from where it follows that the solutions ao lP +i; Pl fe(r) ) ffli p+i ; p fc(7"), a.2,p+i;p,fc(T) are polynomial as 
no (1 ± t) _1 terms have to be integrated. Note also, that the multiplication of J X p ^ p B ' p+ i- p &r 
by X p+ i- p cancels out all rational expressions. I 



7 The order p + 2 



It is now shown how the T 2 k -i+j harmonics at order 0(p p ) in the initial data can give rise to 
solutions to the transport equations with logarithmic divergences at the critical sets 1 ± . The 
main result is the following: 

Proposition 5. There are initial data satisfying the regularity condition \50\) such that the func- 
tions ao >p +2;p,k( T )> a i.p+2-,p,k(T~) and a2, P +2;p,k(i~) with k — 0, ...,2p, solving equations S3 3 a)) , 
S33b\) and S3 3c}) contain logarithmic divergences at t = ±1. 

Proof. A look at the integral formula (|51[) reveals that the order p + 1 solutions ao, p +i ; p,fc(T), 
®i,p+i;p,k(T~) and ci2,p+i;p,fe(T) are composed by two independent parts: one having to do with the 
homogeneous part of the solution which can be eliminated by setting ao,p+i;p,fe(0) — &2,p+i;p,fc (0), 
k = 0, . . . , 2p in the initial data, and the other the nonhomogeneous solution proper which can 
be "turned off" by setting a Oj) . p k(0) = — assuming the regularity condition (j50)) holds. For the 
purpose of exhibiting the existence of logarithmic divergences in the solutions to the order p + 2 
Maxwell transport equations it will suffice to analyse what happens when only the homogeneous 
terms of the order p + 1 solutions are considered — exploiting the linearity of the transport 
equations. Let Y p+l . p k denote these terms. From formula (|51|) one has that 

p h +1 (r) = X p+1 (r)X p ^(0) , (67) 

where again here and in the sequel the subindex string - p ^ has been suppressed. Consequently, 
performing the replacement p — > p + 2 and q — > p in the equations p6ap - (|36b[) and setting 
Q-o,p( T ) = Q-i,p( T ) = o>2,p(t~) = — which is the case if ao, p (0) = — one finds that the components 
bo,p+2 and &2,p+2 of the column vector B p+ 2 containing the nonhomogeneous part of the system 
of ordinary differential equations to be solved are given by 

&o,p+2 = rfq~)3 (^oao.p+i + 52a 2 ,p+i) , (68a) 

&2,p+2 = ^ (ffo a O,p+l + 52 a 2,p+l) , (68b) 

with the polynomials go and 52 given by 

5 o(r) = ^(1 + r) ^(4p 2 + 9 P )t 5 + (8p 2 + Wp - 12)r 4 - {16p 2 + 33p - 4)r 3 

,2 1 iio„ 1 tA_2 //id„2 



Y. 



(56// + U2p + 8)t z - (48// + 126/> + 72)t - (24p + 48) (69a) 



.92 (r) = -^-(1 + r) ( (4p + 9)r 4 - (4p + 6)r 3 



-(28/> + 57)r J + (16p + 24)t + (24p + 42) . (69b) 
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In the above expressions the order p + 1 Maxwell transport equations have been used to eliminate 
the derivatives a' p+1 , s'j p+1 and a 2 p+i that arise. Further, equation (|35p has also been used to 
dispose of ai, p +i. 

For conciseness of the presentation, let 



x- + 1 1 (o)y p+1 (o)= r°). (70) 



Hence, one has that 



Y h _ { a , P+l {r)\_ (zo(l- r)P+V 1 (r) + (-1)p +1 z 2 (1 + r)PJ 3 (r) 1 



a 2 , p+1 (r); V + t) p+2 J?(t) + (-l)^o(l - r)f J|(r) 

where Ji and J3 are the polynomials (|60a |) — (|60bj ) arising from the order p+1 fundamental matrix. 

As in the proof of proposition [4J the crucial point in the argumentation is to analyse the 
decomposition in partial fractions of the entries of 

1 1 [ (1 - r) p+1 ^° a °'P +1 + s ^-p+^ K i + Jl ^ ly+i ( s o q o,p+i + 4 a 2, P +i)K 3 ' 

X p+2,p B P+2 = ~jFd 1 f-l) P 

p+2 > p \ ( 1 + T }P+l ( s o a o,P+i + s s 2 a 2 , P +i)Ki + - Y T y+- A (aoQo.p+i + s 2 a 2 ,p+i)K^ 

(72)' 

where it has been used that go = (1 + t)sq and g 2 = (1 + t)s 2 , and 

^i(r) = ^3^t 3 '- p " 1) (r), (73a) 

^3(r) = ^(l + r^P^-^W, (73b) 

are the Jacobi polynomials arising in the fundamental matrix X p+2jP . 

Contrary to what happens in the analysis of the order the functions ao. p +i and a 2lP +i 

are not homogeneous in (1 ± r) m for some integer m. This fact complicates the analysis of 
the decomposition in partial fractions of the expression (|72p . However, if one is interested in 
exhibiting the generic existence of logarithms in the solutions, their appearance for a restricted 
type of initial data suffices. With this idea in mind, set z 2 = in equation (fTTj) so that 

a , p+1 (r) = *o(l - r) p+2 Jr(r), a 2 , p+1 (r) = z (-l) p (l - rf J|(r). (74) 

So, one finds that in particular 

1 s (-l) p 

-J— ^TT>oao, P+ i + s 2 a 2 ,p+i)-KT + ^ _ r ^ p+3 (so q o,p+i + s^p+i)^ 

( Tj 

where 

G = (-l) p (l - T?Kma 2 + (-1)"(1 - t) 2 K 3 J iS s + K 3 J 3 s s 2 . (76) 

As in the previous section, the problem reduces to calculating G" and finding out whether it has 
zeros at r = 1. 

A direct calculation shows that 

G(l) = p(p + 2) 2 2- 2p - 1 m^+ i 1 ), (77a) 



G'(l) = --(p + 2)(3p J - 4p 2 - 16p + 2) (^ + i j , (77b) 
G"(l) = ^(p + 2)(3p 4 - 4p 3 - 3p 2 - 36p - 98) + . (77c) 
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It can be seen that the only real positive root of the polynomial in p, 3p 4 — 4p 3 — 3p 2 — 36p — 98, 
lies in the interval (3,4) so that one can conclude that G"(l) ^ 0, for p > 1. Using arguments 
similar to those of the proof of proposition 0] one concludes that 

G{t) = (1 - T) 3 h 3 (r) + (1 - r) 2 h 2 (l) + (1 - r)/n(l) + G(l), (78) 

for some polynomials hi, h% and /13, with hi(X) 7^ 0. Whence 

G(r) ~ fta(l) Ml) G(l) , . 

_^_/ l 3( r ) + ^ + _^ + ___. (79) 

Thus, it follows that the integration in formula (|5ip renders logarithmics terms ln(l ± r). Note, 
again, that the multiplication of the integral by X p+ 2,p cancels out all fractional terms in (1 ± r) 
so that the solution lies in the ring generated by polynomials in r and ln(l ± r). ■ 



8 Conclusions 

The present work has been a first attempt to understand, in a general way, the mechanisms 
logarithmic divergences arising at the critical sets I in the solutions to the transport equations 
at X implied by a Maxwell field propagating on a Schwarzschild background. The ultimate goal 
behind this programme is to develop ideas and methods which could be employed in the analysis 
of the much more challenging case of the behaviour of the transport equations implied by the 
Conformal Einstein equations. Preliminary calculations show that by the methods employed in 
this article it should be possible to at least obtain an analogue of propositions [3] and [4] for this 
case. That is, starting from initial data which are exactly Schwarzschildean up to order p — 1, 
and with a first deviation from Schwarzschild data occurring at order p, it should be possible 
to fully integrate the order p transport equations — note that in this case the system consists of 
50 equations! — to show that if a certain regularity condition in the data holds — the analogous 
of condition (|50p — then there are no logarithms in the solutions of the transport equations at 
order p and p + 1. Quite remarkably, some explicit calculations for concrete values of p show 
that there should not be any logarithms in the solution at order p + 2, and that the analogue of 
proposition [5] occurs in this case at order p + 3. It is quite possible, however, that a proof of these 
last two observations is beyond the possibilities of the present methods, and that further ideas 
and insights have to be developed. 

In the light of the evidence gathered in this article one can not avoid to feel that there is some 
powerful group theoretical structure lurking behind which is the responsible of the remarkable 
patterns that have been observed in both the case analysed in this article and more generally in 
the conformal Einstein field equations. 

It would be of particular interest to reformulate the results obtained here in a context or 
language more in tone with ideas and concepts of Differential Galois theory, which in turn would 
render further insights — see e.g. [6] for a brief introduction to the ideas and goals of Differential 
Galois theory. For example, propositions El H] and [5] state that the Galois groups of the order p 
and p + 1 transport equations are trivial, while that of the order p + 2 equations are not. Now, 
the question is the following: is it possible to obtain these results by — say- looking directly at 
the field equations and without having to solve almost explicitly the equations? 
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A Some spinors 



Let {oa, la} denote normalised spinor dyads, cabo l = 1, where cab is the standard alternating 
spinor. In the text, the following spinors have been used: 

taa' = oaoa' + ^A 1 * A' j (80a 

(80b 

z A b = —^o A o B , (80c 



e AB = °AOB, ^AB~°(A^B), ^AB — L A^B, 

1 / \ 1 

XAB = —^(OAtB + tAOB), VAB — ~—j=iA^B, 



e ABCD = OAOBOCOD, ^ABCD 
h-ABCD — —£a(c£d)b- 



l(A°BOCt'D), 



V2 

£ ABCD = L AtBtCt-D, (80d 

(80e 



B The expansions of the Schwarzschild spacetime 

The following expansions of the field quantities on the F-gauge for the conformal Schwarzschild 
solution have been used in the discussion of the article. The expansions have been calculated 
using the computer algebra methods developed in [13] . 

c° = -t + f ~mr 3 - \raT b ^j p + i ^mV + ~mV - 3mV - 2m 2 r 3 ^ p 2 + 0(p 3 ), (81a) 

c 1 = p + i (-4tot 2 + ^r 4 ) p 2 + | (l2m 2 r 2 + 15m 2 r 4 - y mV + ^m 2 r 8 ) p 3 + 0(p 4 ), (81b) 

c± = 1 + f mr 2 - -mr 4 ) p + 1 f — m 2 r 8 - -m 2 r 6 + 3m 2 r 4 - 2m 2 r 2 ) p 2 + 0(p 3 ), (81c) 
Y 6 / 2! y 14 9 / 

/ = 1 + imr 4 p + ^ Qm 2 r 8 - ^m 2 r 6 + ^mV - 2m 2 r 2 ^ p 2 + C(p 3 ), (81d) 
j= = Qr 2 m - imr 4 ) p + | ( - ^mV - ^r 6 m 2 + ^mV) p 2 + G(p 3 ), (81e) 
X 2 = (4mr 3 - 12rm) p + ^ ^24m 2 r - 8m 2 r 3 + 4r 5 m 2 - ^j"™ 2 ^ /° 2 + 0(p 3 ), (81f) 

Xh = h (~ it 3 ™ 2 + r 5 ™ 2 - § mV ) p2 + ° (p3) - (81g) 

For completeness the components of the curvature are also given: 



_± = 4mrp + _ (jt'W _ _ T J m J _ 12m V j ^ + (82a) 

6 2 = 6m (1 - r 2 ) p + ^ ^-12m 2 + 12m 2 r 2 - 10r 4 m 2 + ^™ 2 t 6 ) P 2 + C(p 2 ), (82b) 



(82c) 



e,^^(4r 2 - V + V)p 2 + 0(p 3 ), 

= -6m + (-18m 2 r 2 + 3r 4 m 2 ) p + i ^yT 8 m 3 + 28to 3 t 6 - 90to 3 t 4 + 36to 3 t 2 ^ p 2 + C(p 3 ). (82d) 
Finally, one has that 

/2 — - 

2r 0011 - 4r 1010 - 2£ + - 2V2 X h - V2/, (83a) 

/2 

2(r 110 o+r 0011 ) = ~^- X 2 - 2V2 Xh , (83b) 

/2 

2r 1100 - 4r 0101 - -2f + - 2V2 Xh + V2f, (83c) 

D 

2(r on - r UO o) = 3£- j^X2 - -^-Xh- (83d) 
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C Some selected properties of the Jacobi polynomials 

The following properties of the Jacobi polynomials have been used in the calculations described 
in this article. A good general reference to this matter is the monograph [TO] , 
The Jacobi polynomials are defined by 

P (a,ffl M = n^ + n+l) vW »\ r(a + /3 + n + m+l) /r-lV" 
n lJ nirfa + ^ + n+l)^ U T(a + m + 1) I 2 J ' 1 J 

In the above expression the combinatorial coefficients are understood in the generalised sense 
involving the Gamma function. That is, 



m\ T(m + 1) 



nj r(n + l)r(m — n + 1) 
Derivatives of the above can be rewritten as lower degree Jacobi polynomials via 



(85) 



The polynomials satisfy the symmetry relation 

p(^)(-r) = (-irp(^(r). (87) 
The value of the polynomials at r = ±1 is given by 

^(i)=( n ; a ), ^ ) (-i)=(-i)"( n + /3 )- (88) 
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